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The thermoelectric properties of in plane heterostructures made of Graphene and hexagonal Boron 
Nitride (BN) have been investigated by means of atomistic simulation. The heterostructures consist 
in armchair graphene nanoribbons to the sides of which BN flakes are periodically attached. This 
arrangement generates a strong mismatch of phonon modes between the different sections of the 
ribbons, which leads to a very small phonon conductance, while the electron transmission is weakly 
affected. In combination with the large Seebeck coefficient resulting from the BN-induced bandgap 
opening or broadening, it is shown that large thermoelectric figure of merit ZT > 0.8 can be reached 
in perfect structures at relatively low Fermi energy, depending on the graphene nanoribbon width. 

The high value ZT = 1.48 may even be achieved by introducing appropriately vacancies in the 
channel, as a consequence of further degradation of the phonon conductance. 


I. INTRODUCTION 

Due to its exceptional electronic, optical, thermal and 
mechanical properties, graphene should not only replace 
conventional materials in some existing applications, but, 
above all the combination of all these unique properties 
is expected to inspire fully new applications. [1, 2] How¬ 
ever, one major drawback of graphene for electronic ap¬ 
plications, at least for its use as a switching device, is 
its gapless character [3] that is responsible, e.g., for the 
low on/off current ratio in graphene transistors. [4, 5] It is 
also an obstacle for thermoelectric applications because it 
makes it difficult to separate the opposite contributions of 
electron and hole states to the Seebeck coefficient, which 
is smaller than 100 mV/K in pristine graphene. [6] Ad¬ 
ditionally, the thermal conductivity in graphene is very 
high, even higher than 4000 W/mK for a single layer, [7] 
which is a strong limitation to achieve high thermoelec¬ 
tric figure of merit ZT. Indeed, for small ballistic con¬ 
ductors this quantity is conveniently defined in terms of 
electronic and thermal conductance as 


where P = G e S 2 is the power factor, K = K e + K p 
is the total thermal conductance, T is the absolute tem¬ 
perature, while G e , S, K e and K p are the electrical con¬ 
ductance, the Seebeck coefficient, the electron thermal 
conductance and the phonon thermal conductance, re¬ 
spectively. 

Actually, the thermoelectric figure of merit in graphene 
is limited to ZT < 0.01, leading to poor thermoelectric 
efficiency. However, in spite of these intrinsic obstacles, 
many suggestions of bandgap nanostructuring have been 
taken up, at least at the theoretical level, to enhance the 
thermoelectric performance of graphene and to achieve 
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FIG. 1. Schematic view of graphene/BN heterostructure with 
sub-sections made of BN flakes attached to the main AGNR 
channel. The two leads are made of BN/G/BN ribbons as the 
hybrid parts of the active device. 

ZT values higher than 1. [8] Indeed, since the pioneer¬ 
ing works of Hicks and Dresselhaus, nanostructuring ma¬ 
terials in low-dimensional systems is strongly expected 
to provide higher ZT and thermoelectric efficiency than 
bulk materials, [9] thanks to (i) the enhancement of See¬ 
beck coefficient resulting from size quantization and (ii) 
the reduction of thermal conductance due to interface 
effects. It has been confirmed experimentally in many 
systems. [10] In the case of graphene, the primary expec¬ 
tation of nanostructuring is the bandgap opening to en¬ 
hance the Seebeck coefficient. The reduction of thermal 
conductance is strongly expected too, if it is not paid by 
a stronger reduction of electronic conductance, since all 
parameters entering the expression of ZT are mutually 
coupled and difficult to control independently. 

As the simplest form of graphene nanostructuring, 
graphene nanoribbons (GNRs) were first investigated. 
With finite-bandgap armchair-edges GNRs (AGNRs), 
significant improvement of thermoelectric properties has 
been predicted compared to 2D graphene. [11, 12] How¬ 
ever, in the best case ZT does not exceed 0.35 for the 
narrowest AGNR with only Mcc = 3 dimer lines in the 
width. To enhance further ZT, it is necessary to design 
more sophisticated ribbons likely to degrade strongly the 
phonon contribution to the thermal conductance while 
retaining high electron conduction properties. Several 
works have suggested to use multi-junction GNRs with 
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alternate sections of different width, different chirality, or 
different chirality and width. [12-14] Values of ZT close 
to one have been calculated in such structures. Promis¬ 
ing results were also obtained in kinked (or chevron-type) 
GNRs[15] or even by combining chevron-type geometry 
with isotope cluster engineering. [16] In zigzag-edge GNRs 
of several micron length, edge disorder and extended 
line defects have been shown to enhance ZT to values 
greater than 2 if electron-phonon scattering may remain 
negligible, [17, 18] It has been also predicted that ZT can 
be improved by using graphene nanomesh (GNM), i.e. a 
graphene layer in which a periodic array of nanopores is 
perforated[19] or by designing vertical junctions in mul¬ 
tilayer graphene. [20] 

Recently, a novel form of hybrid monolayer struc¬ 
tures consisting of a mixture of graphene and hexagonal- 
Boron Nitride nanodomains has been proposed and syn¬ 
thesized successfully, [21-23] thus opening new strategies 
for bandgap engineering and designing devices. Some 
theoretical works have demonstrated that Graphene/BN 
heterostructures made of adjacent armchair ribbons are 
always semiconducting, [24-29] while zigzag hybrid rib¬ 
bons may be either semimetallic or semiconducting, 
depending on the type of bonding at graphene/BN 
interfaces, [26, 29, 30] which is promising for transistor 
application with efficient on/off current switching. [31, 32] 
Some works also considered the thermal transport in 
these types of heterostructures [33-35] and have shown 
that they offer the possibility to tune and reduce strongly 
the phonon thermal conductance. Hence, thanks to the 
large bandgap opened in Graphene/BN hybrid structures 
that should lead to good power factor, it is expected that 
high ZT can be achieved. However, according to our best 
knowledge only two works have been reported so far on 
the thermoelectric properties of monolayer graphene/BN 
structures, by Yang et al.[36], and more recently by 
Vishkayi et al. [37] , respectively. The latter focused on hy¬ 
brid zigzag ribbons with a graphene/BN junction and two 
leads made of square lattices. Though an enhancement of 
ZT was observed compared to pure zigzag graphene its 
maximum value was still limited. In the former work, [36] 
the authors reported a significant improvement of ZT up 
to 0.7 for some configurations of superlattice armchair 
ribbons consisting of alternating graphene and BN sec¬ 
tions. However, the insulating effect of BN sections in¬ 
duces a high electron scattering, which shifts the good 
ZT values to high chemical potentials of about 2 eV, i.e. 
a range of energy difficult to exploit in practice. 

In this article, by means of atomistic calculation 
we investigate the thermoelectric performance of GNRs 
with parallel graphene/BN interfaces appropriately dis¬ 
tributed. We demonstrate that by graphene/BN inter¬ 
face engineering it is possible to enhance the phonon scat¬ 
tering in armchair GNRs and thus to reduce strongly the 
phonon thermal conductance while the electronic con¬ 
ductance is weakly affected. It leads to high values of 
ZT at rather small chemical potential. Additionally, we 
show that the figure of merit may be enhanced up to 


1.48 at room temperature by introducing vacancies in 
the graphene region. 

The paper is organized as follows. The model and 
methodology are presented in Sec. II, while the results 
are discussed in Sec. Ill, where we emphasize the role of 
interface phonon scattering (sub-Sec. Ill A) and that of 
vacancies (sub-Sec. IIIB). The conclusion is in Sec. IV. 

II. STUDIED DEVICE AND METHODOLGIES 
A. Device structure 

The structure investigated in this work is schema¬ 
tized in Fig 1. It consists in a graphene ribbon 
to the sides of which BN flakes are periodically at¬ 
tached. It is thus an alternating arrangement of 
graphene and BN/graphene/BN (BN/G/BN) sections. 
Such BN/G/BN ribbons are known to open or broaden 
the bandgap of free GNRs, [25, 27] and are expected 
to modify strongly the phonon dispersion, with thus a 
strong impact on the phonon conductance. In contrast, 
the electronic conductance in the central graphene rib¬ 
bon is expected to be weakly affected by the presence 
of BN flakes. We have focused our study on armchair 
structures that usually provide higher ZT than zigzag 
edge ribbons. [36] 

The width of each sub-region is characterized by the 
numbers of dimer lines Mqc and Mbn- Along the 
transport direction, the structure is characterized by the 
number ubn of BN/G/BN sections and the numbers of 
unit cells N vc and Nbn in graphene and BN/G/BN sec¬ 
tions, respectively. For instance in case of the structure 
sketched in Fig. 1, we have Mcc = Mbn = 4, N vc = 
Nbn = 2, ubn = 2. The total number of unit cells in the 
active region is Na = (ubn + 1) -V„ c + Ubn^bn- This 
active region is connected to BN/G/BN leads having a 
large bandgap. It is assumed that the system can be ex¬ 
tended periodically along the y direction, which requires 
Mcc + Mbn to be an even number. [25] 

B. Methodologies 

The electronic properties were investigated using a 
nearest-neighbor tight binding (TB) approach. The 
corresponding atomistic Hamiltonian takes the general 
form [38] 

*. = £* i*) (*i - E ^ w 01 ( 2 ) 

i (hi) 

where e, is the on-site energy at site i and t, i; - is the 
hoping energy between atoms at i-tli and j-th sites. 
Starting from the second-nearest-neighbor TB model of 
Ref. [25] , we made some small changes to work with a sim¬ 
ple nearest-neighbor model that reproduces very well the 
energy band structure of armchair graphene/BN ribbons 
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FIG. 2. (a) Phonon conductance as a function of temperature for Mcc = 5 (solid lines) and Mcc = 6 (dashed lines), for 
different values of the number nBN of BN/G/BN sections (b) Phonon conductance at room temperature for Mcc = 5 and 
Mcc = 6 as a function of ubn (red squares) and Mbn- (blue circles), (c) Phonon density of states in AGNR with Mcc — 5 
and in BN/G/BN ribbon with Mcc = 5 and Mbn = 9. 


computed using first principle method. [25] The values of 
the parameters used are listed in table I. 

The phonon properties were computed using the fourth 
nearest-neighbor force constant (FC) model. [39] The gen¬ 
eral form of the Hamiltonian for the vibrations of atoms 
was given in [40] and can be rewritten as 

Hp ^ ] —Mi u i iii T ^ ] — (Ui Uj ) (ttj Uj ) (3) 

* <*b> 

where Mi and Ui are the mass and the time-dependent 
displacement of the *-th atom, respectively, and is 
the 3x3 coupling tensor between the i-th and j- th atoms. 
Accordingly, in the small displacement limit, the motion 
equation of the i-tli atom writes[12, 39] 


Mi 


d 2 Ui 

dt 2 


( u i “ u i) 


( 4 ) 


For a finite number of atoms, the system of equations 
4 may be written in a matrix form. After Fourier trans¬ 
form, it can be rewritten as a function of the angular fre¬ 
quency u> of the atomic vibrations as u> 2 U = DU , where 
U is the column vector formed by the relative displace¬ 
ment vectors Ui , and D is the dynamic matrix defined 
as 


jj’i 

-^3x3 

= 

[| 

f - / K ' J for j/* 

\ Mi Mj J 7 

£ Aff for 3 = i 




K. n^i J 


TABLE I. Tight binding parameters for electron study. 


E Ca 

E Cb 

E b 

En 

tcc 

tBN 

tsc 

tNC 

eV 

eV 

eV 

eV 

eV 

eV 

eV 

eV 

0.015 

-0.015 

1.95 

-1.95 

2.5 

2.9 

2.0 

2.0 


It should be noted that the term Y is a 3 x 3 matrix 

n^i 

that plays the same role as the on-site energy in the case 
of electrons. However, this term is varying for each atom, 
depending on the number and nature of couplings with 
other atoms. It is in contrast with the case of electrons 
wherein the on-site energy depends only on the atom on 
this site. The coupling tensor K l:j between the i-th and 
j- th atoms is defined by using a unitary rotation in the 
plane of the mono-layer graphene/BN structure 

Kij = U~ x (dij) K°ijU (Oij) (6) 

where U (Oij) is the rotation matrix[39] and Oij is the 
anticlockwise rotating angle formed between the positive 
direction of the cc-axis and the vector joining the i-tli to 
the j- th atoms. Finally, K°ij is the force constant tensor 
given by 

/ 0 0 \ 

K\j = 0 $ ti 0 (7) 

V 0 0 / 

where <f) r ,$ ti and 4> 4o are the force constant coupling 
parameters in the radial, in-plane and out-of-plane direc¬ 
tions, respectively, depending on the level of neighboring. 
Within a fourth nearest-neighbor, we thus need twelve 
parameters for graphene and fifteen parameters for BN. 
In this work the parameters were taken from Ref. [41] for 
graphene and from Ref. [42] for BN. For the coupling of 
C atoms with B and N atoms we have taken the average 
value of the force constant parameters between graphene 
and h-BN.[36] The masses of C, B, N atoms were taken to 
be equal to 1.994 x 1CT 26 ,1.795 x 1CT 26 and 2.325 x 10" 26 
kg, respectively. 

Due to the similarity between the TB Schrodinger 
equation for electrons and the FC dynamic equation for 
phonons, they can formally be solved using the same ap¬ 
proach. The Green’s function method has been used 
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FIG. 3. (a) and (b) Electrical conductance and (c) (d) Seebeck coefficient for the structures with (a) (c) Mcc = 5 and (b) (d) 
Mcc = 6 for different numbers nBN of BN/G/BN sections. Other parameters: T = 300K. 


in the ballistic approximation. [35] In the case of elec¬ 
trons the total Hamiltonian can be separated in the de¬ 
vice Hamiltonian H R of the active region, the Hamil¬ 
tonian H L ( R ) of the left (right) lead and the coupling 
term t R ( R ) between the device and the left (rigt) lead. 
First, the surface Green’s functions G° L ^ R \ of the un¬ 
coupled leads are calculated using the Sancho’s itera¬ 
tive scheme. [43] They are used to calculate the self¬ 
energies = t l(R)G 0 li ^ r ^ l ^ describing the device¬ 

lead couplings that take place into the device Green’s 
function G = [E + irj — H R — T, S L — £Jj] -1 . In the case 
of phonons a similar expression of the Green’s function 
is found by replacing the energy E of electrons with the 
square of the phonon angular frequency u> 2 . A recur¬ 
sive technique is applied to reduce the size of the device 
Green’s function to the size of the Green’s function Gn 
of the first layer. [44, 45] Then, the electron and phonon 
transmissions, T e ( E ) and T p (w), can be computed as 


T = Tr {ri \i (Gn - Gu*) - GnT^Gn^j } (8) 

where TJ^ = i i s the surface injec¬ 

tion rate at the left (right) lead. Finally, once the trans¬ 


mission T e ( E ) and T p (w) are calculated, the electrical 
conductance, the Seebeck coefficient, and the electron 
and phonon contributions to the thermal conductance 
can be computed as [8] 




Ge (ji,T) 
s (m, T) = 
K e (n,T) 


-- e 2 Lq (h,T) 
Zb (n,T) 
eT Lq(h,T) 


]_ 

T 


L 2 (h,T) 


^i(mT ) 2 

L 0 (n,T) 


I< P = 


r dm 

J 2 ?r 


Huj T p (w) 


dn(cj,T) 

dT 


(9) 


where the intermediate function L n (/q T) is dehned as 


+oo 

r ( 2 [ jvrr rv 1" -9fe{E,n,T) 

Ln Ub T) = - J dET e (E) (E - n) - — - 

— OO 

(10) 

In these equations, f e ( E , /q T) is the Fermi distribution 
function, n(uqT) is the Bose-Einstein distribution func¬ 
tion, T is the temperature and [i is the electron chemical 
potential. More practical forms of the intermediate func- 
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tion and the phonon conductance can be derived as (a) 

{ +oo 

L n {g,T) = { f dE T e (E) (2k b T) n - 1 g e n (E,g,T) 

— OO 

oo 

Kp = ^ J duT p (cj) g p (w, T) 
o 

(U) 

where g e n {E,g,T) = /cosh 2 ( A-=ii), 

g p {oj,T) = ^ 2 T^r) /^inh 2 (/ T^r) and > is the 

Boltzmann constant. Once G e , S , K e and K p are 
obtained from equations 9, the figure of merit ZT can 
be computed from equation 1 


III. RESULTS AND DISCUSSION 
A. Phonon scattering: a key to enhance ZT 

In the following analysis of the thermoelectric perfor¬ 
mance of the device schematized in Fig. 1, unless other¬ 
wise stated we will focus mainly on two cases of graphene 
ribbon width, i.e. Mcc = 5 and Mcc = 6, that corre¬ 
spond to a metallic and a semiconducting ribbon, re¬ 
spectively, in the case of pure AGNR. Similarly, the BN 
width Mbn will be chosen equal to 9 and 12 for Mcc = 
5 and 6, respectively, that correspond to the width be¬ 
yond which the transport properties are unchanged, as 
will be discussed later. The length of active region can 
be adjusted by changing the number ubn of BN/G/BN 
sections, which will be used as a parameter in what fol¬ 
lows. Indeed, by increasing this number one may expect 
to reduce the phonon conductance as a consequence of 
enhanced phonon scattering. 

In Fig. 2(a) we plot the phonon conductance for Mcc 
= 5 and Mcc = 6 as a function of the temperature, 
for different numbers ubn of BN/G/BN sections in the 
active device. We also plot in black lines for compari¬ 
son the phonon conductance of the corresponding pristine 
AGNRs. For both widths we obtain a strong reduction 
of K p in hybrid structures compared to graphene. Ad¬ 
ditionally, it can be seen that the phonon conductance 
decreases further when hbn increases. For instance, in 
the case of Mcc = 5 at room temperature (Fig. 2(b), 
red squares; solid line) i.e. it drops from 0.244 nW/K 
for ubn = 1 to 0.196 nW/K for ubn = 5. Then the 
phonon conductance saturates when the active region is 
long enough because the system tends to behave as an in¬ 
finite superlattice system. Finally, for ubn — 10 at room 
temperature, K p is about 75% smaller than in the corre¬ 
sponding pristine graphene ribbon. Similarly, for Mcc 
= 6 the phonon conductance falls from 0.96 nW/K (pris¬ 
tine AGNR) to 0.346 nW/K (ubn = 1) and even 0.24 
nW/K (ubn = 10). A similar behavior is observed when 
the width Mbn (blue circles in Fig. 2(b) of the BN flakes 
increases. A minimum value of about 10 is required to 
reach the minimum of conductance. 
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FIG. 4. (a)(b) maximum of thermoelectric figure of merit 

as a function of ubn for (a) Mcc = 5 and (b) Mcc = 6. 
Thermoelectric hgure of merit ZT as a function of chemical 
potential g for (c) Mcc = 5 and (d) Mcc = 6. T = 300K. 


This dramatic reduction of phonon conductance is 
not only due to phonon scattering at graphene/BN 
boundaries, [46] but above all to the scattering induced 
by the mismatch of modes between pristine graphene 
and BN/G/BN sections that reduces strongly the phonon 
transmission. This mismatch of modes is illustrated in 
Fig. 2 (c) where we plot the phonon density of states 
calculated in both sections assumed to be infinite. 

The electrical conductance G e and Seebeck coefficient 
S are plotted in Fig. 3 as a function of chemical energy 
g in both cases of Mcc = 5 and 6 and, for compari¬ 
son, in the case of pristine GNR (black dashed lines). In 
Figs. 3(a) and 3(c), i.e. for Mcc = 5 corresponding to a 
metallic AGNR, we can observe the opening of a signif¬ 
icant conduction energy gap, which manifests itself not 
only in the vanishing of G e at low energy, but also in the 
strong enhancement of S, compared to its pristine GNR 
counterpart, whatever the value of ubn- In the case of 
Mcc = 6, for which the pristine GNR is semiconduct¬ 
ing, the bandgap is slightly broadened, which yields a 
small enhancement of the Seebeck coefficient, indepen¬ 
dently of the number of periods ubn (Figs. 3(b) and 
3(d)). It is remarkable that in contrast with the phonon 
conductance, the electronic conductance is not severely 
degraded with respect to the pristine GNR, whatever the 
value of ubn, at least at low energy (first step of trans¬ 
mission) at which most of the conduction takes place. 
Regarding thermoelectric performance, it is definitely a 
strong advantage of this type of nanostructuring design. 
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It should be noted also that the maximum Seebeck coef¬ 
ficient is higher in the case Mcc = 6, corresponding to a 
larger bandgap. This behavior is in agreement with the 
general result that for any material the maximum value 
of S should be proportional to the bandgap, [47] which 
has been observed in several other cases. [12, 48] 

As a consequence of the strong reduction of phonon 
conductance, together with high power factor, good ther¬ 
moelectric figures of merit are expected in theses hybrid 
structures. Since the phonon conductance appears to 
be strongly dependent on the alternating arrangement of 
graphene and BN/G/BN sections in the active region, the 
results should depend also on the lengths N vc and Nbn 
of these sections. Our investigations led us to the con¬ 
clusion that the maximum value of ZT is reached when 
N vc /N bn ~ 1, which is in agreement with the conclusion 
of Yang et al. [36] in the case of a transverse superlattice 
of graphene and BN sections. We found also that for 
Mcc = 5 and Mcc = 6 the maximum of ZT is obtained 
for N vc = Nbn — 8- However, for AIqc = 3 the best ZT 
is achieved for N vc = 6 and N B n = 4. 

In both cases of Mcc = 5 and Mcc = 6, we plot in Fig. 
4 the maximum value of ZT as a function of the number 
ubn of BN/G/BN sections in the active region (Figs. 
4(a) and 4(b)), and the evolution of ZT as a function of 
the chemical potential in the leads, for different values 
of hbn (Figs. 4(c) and 4(d). In all cases, the case of 
pristine GNR of same width is shown for comparison. 
Due to the metallic behavior of AGNRs with Mcc = 
3 p + 2, the AGNR with Mcc = 5 exhibits a very small 
Seebeck coefficient (see Fig. 3(c)) and a ZT of about 
0.03. Thanks to the bandgap opening and the dramatic 
reduction of phonon conductance, the figure of merit is 
strongly enhanced in the hybrid structure (see Figs. 4(a) 
and 4(c)), reaching 0.6 for ubn = 1 and 0.81 for longer 
devices with ubn > 5 . More interestingly, it should be 
remarked that the peak values of ZT are observed here 
for the quite low chemical energy /z = 0.41 eV, which is 
much smaller and more accessible than in the case of the 
transverse graphene/BN superlattice. [36] 

For Mcc = 6, corresponding to a semiconducting 
AGNR, the peak value of ZT reaches 0.65, which is about 
3.25 times greater than in the AGNR counterpart (Fig. 
4b). Due to the larger bandgap, this peak value is ob¬ 
tained for a chemical potential p, = 0.7 eV, i.e. higher 
than in the case Mcc = 5. In the case of the smallest 
graphene ribbon Mcc = 3 with a large conduction gap, 
a very high ZT can be achieved, as shown in Fig 5. The 
maximum value ZT max = 1.14 is reached for ubn = 10 
and the chemical potential g = 1.15 eV. It is again an 
improvement by a factor of 3.25 compared to the value of 
0.35 in the case of pristine AGNR of same width (black 
dashed line). 

The inset in Fig. 5 shows the evolution of ZT as a func¬ 
tion of the width Mcc of the main GNR channel. In prin¬ 
ciple, the maximum of ZT reduces when increasing Mcc 
due to the increase of phonon conductance. However, the 
behavior is not monotonic and the value of ZT is smaller 



FIG. 5. ZT as a function of chemical potential for the thinnest 
ribbon Mcc = 3 in the cases of pure GNR (dashed line) and 
hybrid heterostructure with Mbn = 11 (solid line), M vc = 
6, Nbn = 4, ubn = 10, T = 300 K. Inset: Dependence of 
ZT m a X on the GNR width Mcc depending on whether Mcc 
can be written in the form 3p, 3p + 1 or 3p + 2, with p is an 
integer number. For the inset, Mbn = 10, 11 for even, odd 
Mcc , respectively and M vc = Nbn = 5, ubn = 10. 


for the group of structures with Alec = 3p +1 compared 
to the groups with Mcc = 3 p and Mcc = 3p + 2. 



M(eV) 


FIG. 6. ZT as a function of chemical potential for different 
temperatures ranging from 100K to 400K for the heterostruc¬ 
ture with Mcc = 5. Here AIbn = 9, and ubn = 10. 

In Fig. 6 we plot ZT as a function of /i for different 
temperatures ranging from 100 K to 400 K in the case 
of the structure with Mcc = 5. Starting from room 
temperature (orange curve) with ZT max = 0.81, when 
the temperature decreases to 200 K and 100K, the peak 
value of ZT increases to ZT max = 0.9 and ZT max = 1.01 
, respectively. When increasing the temperature to 400 
K, the figure of merit reduces to ZT max = 0.74. This 
reduction of ZT max when the temperature increases is 
accompanied by a shift of the peak to the left, i.e. to 
smaller chemical potentials. It is due to the broadening 
of the functions g e n (E, g , T) at increasing temperature, 
which moves the highest peak of Seebeck coefficient to 
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FIG. 7. Effect of a single vacancy on the thermoelectric properties for the structure with Mcc = 5. (a) The seven vacancy 
configurations investigated here, (b) Phonon conductance for the different vacancy configurations. ZT and (inset) thermal 
power for the different configurations of vacancy in (c) graphene and (d) BN. Symbols correspond to results without any 
vacancy. T = 300 K. 


lower values of /i. 

B. Effect of vacancies 

In this next sections we will discuss the effects of va¬ 
cancies on the thermoelectric performance of the hybrid 
structure. We will focus our investigations on the struc¬ 
ture with Mcc = 5 for which we observe high ZT at 
relative low chemical potential. 

It has been shown that edge disorder and inner va¬ 
cancies may have a significant influence on the thermo¬ 
electric performance of GNRs.[17, 49] In particular, de¬ 
pending on the vacancy position the figure of merit ZT 
can be larger than in a perfect GNR. Most of works 
on vacancy effects in graphene were focused on either 
random or periodic arrangement of vacancies. To un¬ 
derstand better the influence of vacancies in our hybrid 
structure, we will investigate the effect of the position of 
a vacancy. We assume that it may be possible to control 
experimentally this position by focused electron beam 
technique. [50, 51] We will first consider different configu¬ 


rations of a single vacancy from the middle of channel to 
the edge, and then the effect of a few vacancies. In our 
structure of Fig. 1, the position of a vacancy is deter¬ 
mined by three indices, i.e. (i) n vac that is the position 
of the unit cell to which belongs the vacancy, running 
from 1 to N a , (ii) L i( j_ vac that is the layer index of a 
layer in the cell, equal to 1 or 2, and (iii) m vac that is 
the position of the vacancy in the layer, ranging from 1 
to M, where M = Mcc + 2 Mbn is the total number 
of dimer lines along the width. For the structure with 
Mcc = 5, Mbn = 9, N vc = Nbn = 8, ubn = 10, we will 
consider specifically seven different configurations of a 
single vacancy at m vac = 12,13,14 (vacancy in graphene) 
and 15, 16, 18, 23 (vacancy in BN) in layer L i( i_ vac = 1 
of the cell n vac = 25 which is the first cell of the second 
BN/G/BN section with respect to the left contact. These 
positions are sketched in Fig. 7(a). 

The phonon conductance is plotted in Fig. 7(b) for 
the different positions of the vacancy in graphene (solid 
lines) or in BN region (dashed lines) and compared with 
the case of no vacancy (circles). We observe a strong 
fall of conductance when the vacancy is located in the 
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FIG. 8. Effect of two and three vacancies in the middle of the 
graphene (CC) ribbon on the figure of merit. Mcc = 5 , T = 
300 K. 


center of the graphene ribbon (red line, conhg-1). It is 
even stronger when the vacancy is on the next lattice site 
(blue line, config-2), with a reduction by almost a factor 
of 2 compared to the no-vacancy configuration. When 
the vacancy is at the edge of the graphene ribbon, the 
reduction of phonon conductance is smaller (green line, 
config-3). When the vacancy is located in the BN region, 
we still observe a reduction of phonon conductance, but 
more limited than when it is in graphene and less sensi¬ 
tive to the position of the vacancy. 

In Fig. 7(c) and 7(d) we plot the thermoelectric figure 
of merit ZT as a function of /i for the different configura¬ 
tions of the vacancy in graphene and in BN, respectively. 
The corresponding power factor P is displayed in inset. 
In comparison with the case without vacancy, we observe 
a remarkable enhancement of ZT up to 1.26 if the va¬ 
cancy is in configuration 1, while we see a degradation 
of ZT in the other two cases of vacancy in graphene. In¬ 
deed, in spite of the strong reduction of phonon conduc¬ 
tance, the configuration 2 surprisingly leads to smaller 
ZT , which is due to a very small power factor, as shown 
in inset of Fig. 7(c). In contrast, the power factor is 
not at all affected when the vacancy is in the center of 
the graphene ribbon, which is also quite surprising. This 
behavior can be understood from the analysis of the den¬ 
sity in the ribbon that is strongly dependent on the lat¬ 
tice site position, as shown in previous works. [32, 52] By 
analyzing the current density in AGNRs Whilhelm et al. 
shown that as a consequence of quantum confinement 
in the transverse flow direction, streamlines appear with 
a threefold periodicity in transverse site position, sepa¬ 
rated by stripes of almost vanishing flow. [52] We observe 
a similar behavior in the graphene part of the structures 
studied in the present work. Starting from an edge of the 
GNR, i.e. m = 1 at the edge, the density is minimum for 
the positions m that are multiple of 3. In the particular 
case of the ribbon of width Mcc = 5 studied here, the 
density in the center of the ribbon (m = 3) is several or¬ 
der of magnitude smaller than at other positions m = 1 
(maximum) and m = 2. Hence, if there is a vacancy at 


this position, i.e. in configuration 1, the effect on the 
density and the conductance of the ribbon is negligible, 
which explains that the power factor is the same as in in 
the no-vacancy case, leading to an enhancement of ZT 
as a consequence of the reduced phonon conductance. In 
contrast, if the vacancy is at any other position in the 
GNR, the density and the conductance are strongly af¬ 
fected, which also degrades ZT in spite of the reduction 
of phonon conductance. 

When the vacancy is in the BN region (Fig. 7(d)) the 
power factor is almost not influenced at low energy, which 
results in a small enhancement of ZT, as a consequence 
of the reduction of phonon conductance observed in Fig. 
7(b). In this case, we always observe an enhancement of 
ZT for any vacancy in a BN region. This is an interest¬ 
ing property of G/BN heterostructures compared to pure 
graphene ribbons in which a vacancy may either degrade 
or improve ZT. 

To fully exploit the influence of vacancies in the center 
of the graphene ribbon, we have made simulations with 
two or three vacancies randomly distributed in the center 
line of the structure with Mcc = 5 and ubn = 10- As 
shown in Fig. 8, ZT can be further enhanced and reaches 
1.48 in the case of three vacancies, as a consequence of 
additional degradation of phonon conductance. 


IV. CONCLUSIONS 

In this article we have demonstrated excellent ther¬ 
moelectric properties of in-plane heterostructures based 
on armchair GNR with appropriate engineering of 
graphene/h-BN interfaces. By attaching periodically BN 
flakes to the GNR, it is indeed possible (i) to reduce 
strongly the phonon conductance compared to pure GNR 
and (ii) to open or broaden the conduction gap which en¬ 
hances the Seebeck coefficient. Additionally, the electron 
conductance is remarkably weakly affected. It finally re¬ 
sults in a large enhancement of the thermoelectric figure 
of merit ZT. In the case of perfect ribbons, the maximum 
value ZTmax of 1.14 is reached for the thinnest ribbon 
with Mcc = 3 carbon dimer lines along the width. For 
larger ribbons with Mcc = 5 and Mcc = 6 we can obtain 
ZT max of 0.81 and 0.65, respectively. Interestingly, in the 
case Mcc = 5 the peak of ZT is reached at relatively low 
chemical potential /x = 0.41 eV. When reducing the tem¬ 
perature from 300 K to 100 K, ZT max increases from 
0.81 to 1.01. The investigation has also shown that engi¬ 
neering vacancies in graphene or BN regions may allow 
to further enhance ZT thanks to additional reduction of 
phonon conductance. However, the electron conductance 
is very sensitive to the vacancy position in graphene. The 
vacancy position must be chosen carefully at a lattice site 
where the electron density is small to avoid unacceptable 
degradation of the power factor. In the case Mcc = 5, it 
is shown that by generating some vacancies in the center 
of the GNR, i.e. where the charge density is negligible 
in a perfect ribbon, the power factor is not affected and 
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ZT max reaches up to 1.48 at room temperature. 
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